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Abstract 

Imagine a swarm of free particles near a point P outside a gravitating mass M and a free reference 
particle at P that is on a radial escape trajectory away from M. Relative to this reference particle 
and in a Fermi normal coordinate system constructed along its worldline, the particles in the 
swarm that move along the radial direction and are ultrarelativistic (that is, they have speeds 
above c/v2) decelerate toward this terminal speed. On the other hand, the swarm particles that 
are ultrarelativistic and move in directions normal to the radial (jet) directions accelerate to almost 
the speed of light by the gravitational tidal force of the mass M. The implications of these effects as 
well as the influence of the higher-order terms on the tidal acceleration mechanism are investigated. 
The observational evidence in support of these general relativistic effects is briefly discussed. 

PACS numbers: 04.20. Cv, 97.60.Lf, 98.58.Fd, 98.70.Sa 
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I. INTRODUCTION 



The geodesic deviation equation in general relativity is often compared with the Lorentz 
force law in electrodynamics. Both connect the corresponding fields to the mechanics of test 
particles and can thus be employed to provide operational definitions of the gravitational 
and electromagnetic fields, respectively. In terms of the dynamics of ultrarelativistic test 
particles, however, there is a fundamental difference between the two equations. To illustrate 
this difference, let us consider the motion of a test particle of mass m and charge q in an 
electromagnetic field (E, B) in Minkowski spacetime with inertial coordinates = (t,x), 

x(^=)=^(E + vxB), (1) 

where c = 1 throughout this paper. This equation can be written as 

^ = — VT^y2 [E - (E • v)v + v x Bl. (2) 

dt m 



It is evident from this relation — as a result of the appearance of \J\ — v 2 = I/7 as a factor 
on its right-hand side — that as v — > 1, it becomes very difficult to change the velocity of the 
test particle; that is, either the field should be maintained over a long integration time or the 
external electromagnetic field must become exceedingly strong. Otherwise, charged particles 
with v 2 extremely close to unity simply behave as essentially free particles and travel along a 
straight line. This may be interpreted in terms of the energy of the ultrarelativistic particle, 
i.e. d(in^)/dt = gE ■ v. For v 2 — > 1, the particle has enormous energy and any significant 
change in its velocity produces a significant change in the energy of the particle; therefore, 
unless an interaction with enormous energy is involved, the particle can be considered to be 
essentially free. Moreover, in the limit of a massless particle such that q/m remains finite 
as m — > 0, the particle simply follows a straight line at the speed of light, as expected. This 
situation should be contrasted with the well-known phenomenon of bending of a light ray in 
a gravitational field; therefore, the analogue of equation (J2J) must be fundamentally different 
for the gravitational interaction. In fact, we will demonstrate that the analogue of the I/7 
factor is absent in the gravitational case. For one-dimensional motion along a straight line, 
(J2J) reduces to 

dv q , 9 . 3 _ 

Tt = ^-v 2 m, (3) 

where E\\ is the component of the electric field along the direction of motion of the particle. 
The relevant factor in this case is I/7 3 ; again, we will show that such a factor is remarkably 



absent in the gravitational case. 

Consider now the corresponding situation in a gravitational field. The analogue of equa- 
tion (JTJ) is the geodesic deviation equation. To express this equation in a form that can 
be compared with (JTJ, we must establish a Fermi coordinate system about one of the 
geodesies jl|. In the curved spacetime of general relativity, such a coordinate patch is the 
closest analogue of the global inertial system employed in equation (JTJ). Let = (T, X) 
be the Fermi coordinate system that is valid in a cylindrical spacetime region around the 
worldline of the reference geodesic observer O. The motion of any other free test particle A 
in this neighborhood is given by the geodesic equation in Fermi coordinates 

d 2 X» „„ dX a dX> 3 



ds 2 + r ^ ds ds - °' (4) 
where s is the proper time of A and T^g is the spacetime connection in Fermi coordinates. It 
is important to recognize that the reference geodesic remains fixed throughout our analysis. 
This approach to the deviation equation should be distinguished from the traditional one 
based on the Jacobi equation, which involves a certain linearization about the reference 
geodesic. Rather, we simply solve the geodesic equation in a Fermi normal coordinate 
system established along the reference geodesic. 

Let us note that equations (HJ-© are expressed with respect to the coordinate time t of 
the global inertial system; therefore, it is necessary, for the purposes of comparison, to write 
equation (jlj) in terms of the temporal Fermi coordinate T as follows: 

fr dx«dx^_ 

ds 2 + r ^ ds ds ~ °' (5) 

r dx 1 dx^_ 

ds 2 + Laf3 ds ds ~ U ' {b) 



Using the identity 



d 2 X i d 2 TdX l fdT\ 2 d 2 X i 

+ hr: -^r> ( 7 ) 



ds 2 ds 2 dT \ds J dT 2 



equations ©-(J7J) imply that 



d 2 X i (■ ,dX i \dX a dX 13 



dT 2 V aP a(3 dT J dT dT 
This is the analogue of equation (JTJ) in the theory of gravitation, except that the geodesic 
worldline could in general be timelike, null or spacelike. To ensure that we are dealing with 
a timelike worldline, equation must be supplemented with the requirement that 

dX^dX u 
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which is preserved by (@J) throughout the motion. Writing the four-velocity of A in Fermi 
coordinates as 

tp = r(i,v) (io) 

with V = dX/dT, equation Q can be expressed as 

T- 2 = -g 0O -2g 0i V i -g ij V i V j . (11) 

Here T = dT/ds is the modified Lorentz factor of particle A, so that r — > oo indicates that 
ds — > 0, i.e. the speed of particle A approaches the local speed of light and the worldline of 
A approaches a null geodesic. 

The Christoffel symbols in equation (jHJ) are obtained from the metric tensor in Fermi 
coordinates 

9oo = -1 - F ll m X i X* + 0(|X| 3 ), (12) 
got = ~ F R 0j i k X j X k + O{\X.\ 3 ), (13) 

9ij = S tj - ~ F R ikj iX k X l + 0(|X| 3 ), (14) 

where 

F Ra/3~fs(T) = _R^p CT A Ai (Q , ) A iy (/3) A p (7) A CT ((5) (15) 

is the projection of the Riemann tensor on the orthonormal parallel-propagated tetrad frame 
A M / <. along the worldline of O. Here A^^ = dx^/dr is the four- velocity vector of O and 
constitutes the temporal axis of its local frame while A M ^, i = 1,2,3, are the unit spatial 
axes that constitute the orthonormal spatial triad of the local frame such that 

g^ (a) y w = r] al3 . (16) 

Here r] a p is the Minkowski metric tensor with signature +2 and r is the proper time of 0\ 
moreover, the Fermi system is so constructed that for this reference observer (T, X) = (r, 0). 

The infinite series in (fT2*j) - (jl4j) consist of increasing powers of the relative distance with 
time-dependent coefficients that involve derivatives of the Riemann tensor along the world- 
line of O. The general behavior of the higher-order terms has been discussed in j^, in 
fact, though the higher-order terms in (|12j ) - (jl4j) can be computed in principle, this has ac- 
tually been done only up to the terms of order |X| 4 . The Fermi coordinates are admissible 
for |X| < TZ, where TZ(T) is the radius of the cylindrical region along the worldline of O. 



To determine TZ, consider the nonzero components of the Riemann tensor and its covariant 
derivatives along the reference geodesic 0\ then, 1Z is defined to be 



inf { 



(3-yS,p 



F R 



...}. 



(17) 



Raf3"/S,p 

The relationship between TZ and the radii of convergence of the series in (fT^ - fpPfj) is not 
known. 

Using the explicit form of the metric tensor in (|12 p -([14 )) . equations (jHJ) and (|TT|) can be 

expressed as 



cPX i 

It 1 



+ F Ro ioj X j + 2 F R ikjQ V k X i 



+ ( 2 ^okjoV'V" + ^ F R ikj iV k V l + 1 F R okj iV i V k V l ) X* + 0(|X| 2 ) = 0, (18) 



and 



J_ = l-V* + F R Qi0j X i Xt + ^ F R Qjik X i Y i X'>- 
1 



+^ F R ikj iV i X k V j X l + 0(|X| 3 ) > 0. (19) 

Equation (jTHjl should be compared and contrasted with equation (J5J): There is no impedi- 
ment in the gravitational case for changing the velocity of an ultrarelativistic particle with 
|V| initially very near unity in contrast with electrodynamics. Moreover, it follows from 
equation (jSJ) that in the gravitational case, the analogue of the electromagnetic particle 
energy equation takes the form 

1 dT _ dX a dX 13 

YdT ~ al3 dT dT ' ( ° } 

so that unlike the electrodynamic case, the rate of variation of particle energy is in fact 
proportional to the energy of the particle. For one-dimensional relative motion in the X 
direction with V = dX/dT, equation (fTHj) reduces to 

^ + K (l_2\/ 2 )X + O(X 2 ) = 0, (21) 

where k = f Rtxtx(T). Equation (|2Tj) should be compared and contrasted with equation J3J). 
Neglecting higher-order terms in ()21|). we note the existence of a critical speed V c = l/v2 ~ 
0.7 in the gravitational case. The physics of this critical speed has been discussed in our 



previous work ; indeed, the existence of the critical speed is basically due to the 

fact that the motion of A is expressed with respect to the Fermi coordinate time T, which 
reduces to the proper time of the reference observer O along its worldline. 

Consider the solution of (|21|) with initial conditions that at T = 0, X = and V = 
^/rjj — 1/To such that To >> 1. Ignoring higher-order tidal terms, a simple integration of 
(|21|) demonstrates that depending on the sign and magnitude of k, one can raise or lower 
the velocity V considerably. In fact, V can be raised to a value such that the quantity T~ 2 
corresponding to equation (|T9j). i.e. 

T~ 2 = l-V 2 + kX 2 + 0(X 3 ), (22) 

vanishes provided that higher-order terms are neglected. This singular behavior signals 
the breakdown of the test particle approach adopted in our treatment, where higher-order 
tidal terms are perforce ignored. These issues will be illustrated in the next section in 
the context of the simplest gravitational field of astrophysical interest, namely, the exterior 
Schwarzschild spacetime. It turns out that the results are observationally significant only 
near gravitationally collapsed configurations. 



II. TIDAL EFFECTS IN SCHWARZSCHILD SPACETIME 

Imagine a swarm of relativistic particles around a point P in the exterior Schwarzschild 
spacetime that corresponds to the gravitational field of a mass M embedded in an otherwise 
flat spacetime with asymptotically inertial coordinates (t, x, y, z). Let the radial line joining 
the center of the spherical symmetry to P be the z axis. Then in terms of the corresponding 
spherical polar coordinates (r,9,<p), the Schwarzschild metric is 

_ rfs 2 = _ (1 _ 2GM } dt2 + (1 _ 2GM^ dr2 + r2( ^ 2 + sin 2 #2) (23) 



It is convenient to refer the motion of the particles in the swarm to an observer O that 
starts from P at r = r > 2GM and moves relatively slowly along the radial z direction 
reaching infinity with a speed that is very small compared to unity. Our results turn out to 
be essentially independent of such a speed; therefore, we set it equal to zero for the sake of 
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FIG. 1: Plot of the initial position and velocity of a particle in the Fermi coordinate system. 
The dashed lines represent the critical velocity cone demarcating the acceleration domains from 
the deceleration domains. The critical angle 9 C is given by tan# c = 1/V C = y/2 , so that 9 C is 
approximately 54.7°. 

simplicity. The geodesic path of O is thus given by 



dt 2GM x dr /2GM 



Tr = {1 -— > ' * = V— ' (24) 
We then set up a Fermi coordinate system (T, X, Y, Z) about the worldline of O such that 
A M ( 3 ) is along the radial direction. In fact, in the (t,r,0,(f)) coordinates we choose 



, I2GM / 2GM\-i 



(25) 



There is then a rotational degeneracy in the choice of A^ and A M , 2 ^; once these vectors are 
chosen at a given instant of time, they are then parallel transported along the worldline. 
The nonzero components of the Riemann tensor along the path of O are given by 

F -Roi01 — F Ro202 = — ^0303 = k, (26) 

F ^2323 = ^3131 = F R\2\2 = ~k, (27) 

except for the symmetries of the Riemann tensor. Here k is given by 

= ^> (28) 



which can be expressed via the integration of (jzljl and r^Tas 

k(T) = ~2GM(rl' 2 + -V2GMT)- 2 . (29) 



2 



for an arbitrary timelike geodesic in 



More generally, the evaluation of equation 
Schwarzschild spacetime has been considered in 

The equations of motion of the particles in the swarm relative to O are thus 

X - \kX[l - 2X 2 + jj(2F 2 - Z 2 )} + ~kX(5YY - 1ZZ) = 0, (30) 

Y - -kY[l - 2Y 2 + -(2X 2 - Z 2 )] + -kY(5XX - 7ZZ) = 0, (31) 
2 3 3 

Z + kZ[l - 2Z 2 + -(X 2 + Y 2 )} + -kZ(XX + YY) = 0, (32) 
3 3 

together with the timelike condition, namely, that 

T- 2 = l-V 2 -h(X 2 + Y 2 -2Z 2 ) + h^XZ-XZ) 2 + (YZ-YZ) 2 -2(XY-XY) 2 ] (33) 

be positive. Here the overdot denotes differentiation with respect to the Fermi time T; 
moreover, we have limited our attention to terms given explicitly in equations (fT2*|) - (fTH) . 
Equations (J30|) - (j32|) and (J33j) can be put in dimensionless form if all spatial and temporal 
durations are expressed in units of GM. There is therefore a simple scaling law at work 
here for different mass scales (for instance, microquasars — > quasars). Our results are most 
important very near the source, since the tidal forces decrease as r~ 3 away from the source. 

To simplify matters, we take advantage of the axial symmetry of this system about the 
Z direction and set Y(T) = for all T. The resulting system can be integrated with the 
initial conditions that atT = 0,X = Z = and 

X = Vocosi, Z = Vosmi, (34) 

where i is the inclination angle and V is the initial speed V = a/Fq — 1 /r such that 
T > 1, see Fig.BQ 

For the motion of a free test particle A along the (radial) Z direction, i = tt/2, the 
equations of motion (|3T)|) - (|3*2*|) reduce to the nonlinear equation 

Z + kZ(l -2Z 2 ) = 0. (35) 

The main aspects of nonlinearity here that are significant for our discussion are the existence 
of critical solutions involving uniform motion at V c = ±1/a/2 and the effective change in 
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FIG. 2: Plot of the Lorentz factor T = l/\Jl- Z 2 - 2GMZ 2 /r 3 versus T/(GM) based on the 
integration of equation (£SJ) for r = 100 GM and T = 10, 20, 50, 100 and 1000. 

the sign of curvature k for Z 2 > 1/2. For initially ultrarelativistic motion (i.e. Vq > l/y/2), 
the particle decelerates toward the terminal speed 1 / y/2 as demonstrated in Fig. |21 The 
details of the deceleration process depend upon r > 2 GM. However, it is important to 
point out a general feature of the tidal deceleration mechanism in the particular case under 
consideration in Fig. [21 the particle starts at T = with ro = 100 GM, Z = and Z > l/y/2 
and decelerates to T < 40 during a time T « 80 GM regardless of the initial r > 50. Thus, 
for r — > oo, "infinite" deceleration can occur in a relatively short period of time; however, 
in this limiting case the test particle approximation breaks down. 

On the other hand, for motion along the X direction, 2 = 0, the equations of motion 
reduce to an equation similar to ()35|). 



X - -kX(l-2X 2 ) 
2 K ' 



0. 



(36) 



except that an initially ultrarelativistic particle accelerates along directions normal to the jet 
direction as illustrated in Fig. El As before, the particle starts at T = with ro = 100 GM, 
X = and X > l/y/2; the integration ends at T end when r — » oo. We note that T end 
may occur outside the domain of validity of Fermi coordinates; however, as To — > oo, 
X en d = X(T en d) ~ Tend tends to zero and one might expect that in this case the influence 
of the higher-order terms may be small. It is in fact possible to show that in this limit T en ^ 
approaches zero as Tq 2 ^ 3 . This is demonstrated in appendix A. It is necessary to remark 
here that for sufficiently large r , the influence of the particle on the background geometry 



can no longer be neglected and hence our treatment becomes invalid. 

To discuss the physics of the tidal acceleration/deceleration phenomena, it is useful to 
have an invariant measure of the energy of the test particle A. Let us therefore introduce a 
class of static observers with four-velocity in the Fermi coordinate system such that in 
the (T, X, Y, Z) coordinates 

£//(T,X)= i- 1 L r , 0,0,0). (37) 

These fundamental observers associated with Fermi coordinates are in general accelerated; 
of course, one exception is the reference geodesic O that is at rest at the spatial origin of 
the Fermi system. At any given event along the path of A, the corresponding fundamental 
observer at that event measures the energy per unit mass of A to be V = —g^U^U 11 , or 

V^o-^LV. (38) 



In the approximation scheme of this section, we have T = ^—goo T, so that with 

-^oo=l + ^(X 2 + F 2 -2Z 2 ), (39) 

T decreases in the deceleration case (X = Y = 0) and increases in the acceleration case 
(Y — Z — 0). Therefore, employing the invariant f instead of T would simply enhance these 
main results of the present section. 

The integration of equations (}3*3*|) in the (X, Z) plane shows that deceleration is maximum 
at i = 7r/2 and monotonically decreases with decreasing inclination until about 35° as 
demonstrated in Fig. HJ Moreover, acceleration is maximum at i = and monotonically 
decreases away from the X demonstrated in Fig. 03 until about 35° when it turns 

into deceleration. This circumstance can be illustrated by means of the critical velocity cone 
in Fig. H a particle with velocity within the cone decelerates relative to O, while a particle 
with velocity outside the cone accelerates relative to O. 

In this section, we have studied the solutions of the geodesic equation limited to the 
lowest-order tidal terms in an appropriately chosen Fermi coordinate patch. Of course, the 
equations of general relativity can be expressed with respect to any admissible system of 
coordinates. We choose a quasi-inertial Fermi normal coordinate system due to its corre- 
spondence with the analysis of observational data. Inside the Fermi coordinate patch, the 
geodesic equation is exact within the test particle approximation scheme. Thus the particle 
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FIG. 3: Plot of the Lorentz factor V = 1/y 1 - X 2 + GMX 2 /r 3 versus T/(GM) based on the 
integration of equation (J3HJ) for r = 100 GM and T = 20, 50, 80, 100, 200 and 300. 



velocity V can take any value consistent with the equation of motion (fTHj) and the time- 
like condition (|19|). see |5[. The Fermi coordinate system has been employed extensively in 
general relativity theory; for instance, it is used in to discuss the local bending of light 
in a gravitational field. The robustness of our numerical results should be emphasized. In 
fact, our results do not change appreciably if instead of X(0) =0, the initial position of the 
test particle is chosen reasonably close to the reference observer. Moreover, the reference 
observer can be any relatively slow-moving test particle on a radial escape trajectory 



III. TIDAL ACCELERATION 



This section is devoted to the singular phenomenon of tidal acceleration of particles to the 
local speed of light within the Fermi coordinate system as shown in Fig. El This singularity 
comes about due to the nonlinear character of equation (JHrJ|) : for V > 1/V2, the nonlinear 
factor 1 — 2V 2 changes sign thereby leading to tidal acceleration. As r increases, T en d/ (GM) 
decreases; the graph of T e „d/(GM) versus log 10 r is given in Fig. El It follows from this 
figure that the singularity under consideration here is dynamic in origin and has nothing to 
do with the kinematic breakdown of the Fermi coordinate system for |X| > 1Z. 

Imagine the general solution of the geodesic equation in the Fermi coordinate system 
in the (X, Z) plane, since the symmetry of the configuration permits us to set Y(T) = 
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FIG. 4: Plot of GM/Tfec versus inclination angle i (measured in degrees in this figure) based 
on integration of equations ()30|) - (|32j) with initial data X = Z = Y = 0, X = Vocosi, Y = 0, 



Z = V smi and r = 100 GM at T = 0, where V = - 1/T with T = 100. The quantity 
Tdec is defined to be the duration of deceleration from Tq = 100 to T = 80. 
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FIG. 5: Plot of GAf/T en d versus inclination angle i (measured in degrees in this figure) based 
on integration of equations (jSQ|) — ({HZ}) with initial data X = Z = Y = 0, X = Vocosi, Y = 0, 
Z = Fosini and r = 100 GM at T = 0, where V = ^T' A - l/r with T = 100. 
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FIG. 6: Plot of T en d/(GM) versus log 10 (ro) based on the integration of equation with ro 



100 GM, X = and X = ^fT 2 - 1 /T Q at T = 0. 



for all T. In this case, the generalization of the system (}3Tjj) - (f3*2*j) that takes into account 
second-order tidal terms |8] is given by 



X + \X(1 - 2X 2 - -Z 2 ) + —XZZ] = \2XZ(1 - 2X 2 - -Z 2 ) 

r 3 L v 3 ; 3 J 2r 4 1 K 2 J 

- (2X 2 -5Z 2 )XZ], 
Z - [Z(l - 2Z 2 + \x 2 ) + \XXZ] = ™ [6X 2 - 12(1 - 2Z 2 )Z 2 



(40) 



- 2QXXZZ - 5Z 2 X 2 - 11X 2 Z 2 ] , (41) 
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where we have employed the linear approximation scheme described in appendix B of p for 
the evaluation of the covariant derivatives of the Riemann tensor along the reference geodesic 
O. In fact, the relevant nonzero covariant derivatives that are needed in the derivation of 
equations (|4(J|) and (|41j) are given by 



^1313,3 — — F -Roi01,3 — — ^0103,1 — -^0303,3 — 3 — , (42) 



except for the symmetries of the curvature tensor. The timelike condition in this case takes 
the form 

-i = l-X 2 -Z 2 + ^\X 2 -2Z 2 --(XZ-ZX) 2 } 
T 2 H L 3 

-^f[(3X 2 -2Z 2 )Z-^(XZ + ZX)XZ 2 ]. (43) 
Moreover, f = y/— g 00 T in this approximation, where 



9oo = l + ™{X 2 - 2Z 2 ) - ^(3X 2 - 2Z 2 )Z. (44) 



We integrate equations (}4T)|) and (|4*T|) with initial conditions that X = Z = 0atT = 
with X = Vq and Z = 0; then, we use the results to plot the behavior of T and f versus 
T/(GM). These plots turn out to be indistinguishable from figure EJ since the contribution 
of the second-order tidal terms is very small; in fact, for r = 50, T en d/ (GM) turns out to be 
72.14 for T and T compared with 72.12 in figure El Thus, the singularity in T is moderated 
by the presence of the higher-order terms. 

These results follow from the integration of equations that are not the actual geodesic 
equation; indeed, the geodesic equation would contain the infinite set of higher-order tidal 
terms, whereas we have considered only the first and the second terms of this infinite set. 
Though the influence of each higher-order term by itself may be small, as we have just 
demonstrated in the case of the second-order terms, we cannot conclude that the same is 
necessarily true of the whole set. Physically, however, our results indicate that with respect 
to the ambient medium around the central source, the free particles of the swarm can undergo 
significant but finite tidal accelerations and decelerations depending upon their directions 
of propagation. The observational aspects of these results are discussed in the following 
section. 
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IV. DISCUSSION 



In our work on the acceleration/deceleration phenomena, we have adopted an approach 
based on relative motion within a quasi-inertial reference frame (i.e. a Fermi normal coordi- 
nate system); in fact, this treatment is generally consistent with the analysis of observational 
data. Moreover, we have concentrated on geodesic, that is, force-free, flow for the swarm of 
particles. We have neglected plasma effects in this paper for the sake of s imp licity; however, 
a more thorough investigation should certainly take these into account jlli I12I ] . Indeed, in 
the comparison of our theoretical results with observation, the astrophysical environment 
around the central source should play an important role. 

It should be emphasized that the acceleration/deceleration phenomena that we have 
discussed occur within the Fermi coordinate system and are thus measurable relative to the 
ambient medium surrounding the central source. Let us note that at T = and r = r , 
where the worldlines of O and A intersect, 

A" (0) ^ = -r , A M (3) f/ M = r Vo sin*. (45) 

These invariants can also be computed in the background Schwarzschild coordinate system 
and we find after a straightforward calculation that 



2GM 

£=(1 + J Vosini)r , £ = r r V r ocosi, (46) 

V ^0 

where £ and £ are respectively the energy and orbital angular momentum per unit mass 
of the particle A as determined by the static inertial observers at infinity in the exterior 
Schwarzschild spacetime. In fact, £ and L are constants along the geodesic orbit of A. Thus 
if this geodesic orbit actually escapes to spatial infinity {£ > 1), it would not appear to have 
any remarkable features as determined by the asymptotically static inertial observers; for 
example, for the accelerating particle discussed in section ITTT1 £ = T and £ = r r Vo since 
i = 0. On the other hand, if the motion of A is referred to the ambient medium surrounding 
the central source, the acceleration/deceleration phenomena that we have discussed can be 
measured. This is the key point: the gain or loss of gravitational tidal energy is measurable 
relative to the ambient medium. In fact, the interaction of (charged) particles in the swarm 
with the ambient medium can lead to the radiation of tidal energy to distant observers. For 
instance, the collisions of such accelerated particles with those of the ambient medium can 
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transfer their tidal energies to the latter particles that could then escape from the system 
and appear far away as highly energetic cosmic rays. 

The acceleration phe nomena are consistent with the recent observations — by the Chandra 
X-Ray Observatory |l3[ — of accelerated motion normal to the jet directions in four neutron 
stars in our galaxy: Crab Pulsar, Vela Pulsar, PSR B1509-58 and SNR G54.1+0.3. A de- 
tailed analysis of the recent Crab nebula X-ray data is contained in 14, 3|. The deceleration 



phenomena are consistent with observations of the speed of the plasma clumps in micro- 

n 

quasar jets [16]. Our theoretical assumptions rely upon the analysis of such observations 
based on the detection of relative motion using the standard flat geometry of an inertial 
system of coordinates. Indeed, the acceleration of particles normal to the jet direction is 
measured relative to the fixed central features associated with the jets near the source jl^ . 
Moreover, the motion of a clump within a jet is measured relative to certain "fixed" features 
of the ambient medium ^J. 

What is the source of energy for the gravitational acceleration of particles? While a local 
description of gravitational energy is not possible in accordance with Einstein's principle 
of equivalence, we must nevertheless account for the energy of the swarm of particles as 
measured by the ambient medium around the source. Ultrarelativistic particles moving 
away from the source within the critical velocity cone lose energy and approach the critical 
speed regardless of their initial ultrarelativistic speed. On the other hand, ultrarelativistic 
particles moving outside the critical velocity cone gain energy and accelerate as they move 
outward. We expect that there is essentially a balance between the energies that are lost 
and gained in the deceleration and acceleration processes, respectively. If the net energy 
that leaves the system is positive, one may speculate that the whole system of particles 
around the source would slightly shrink. Such a contraction would lead to the release of 
gravitational energy that could account for the net loss of energy to the highly energetic 
particles that leave the system; however, the detailed dynamics as well as the limitations of 
such processes is beyond the scope of this work. 

The rotation of the central source has been neglected in this paper; however, detailed 
studies have revealed that — other than specifying the jet directions — the influence of the 
rotation of the central source on our specific acceleration/deceleration results is rather 

DQDDD 

Finally, we should mention that UHE cosmic rays of energy ~ 10 20 eV are not expected 
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Our ac- 



to reach the Earth from extragalactic sources due to the GZK effect \v\ . 18. \v\. 
celeration results indicate that the observed UHE cosmic rays may come from microquasars 
or neutron stars in our galaxy. This idea could be tested with the Pierre Auger Observa- 
tory [20]. 



APPENDIX A 

The purpose of this appendix is to study the behavior of T en d as r — > oo; indeed, we will 
show that in this limit T en d oc T^ 2 ^ 3 . As in the numerical work reported in this paper, we 
will consider the solution X(T, Vo) of the differential equation 

X = - ( =^Lx{l-2X 2 ) (Al) 

with initial conditions X(0, Vo) = and V(0, Vo) = Vo, and 

^=l-f 2 + ^ 2 - (A2) 

We note that the solution of the differential equation has continuous partial derivatives of 
all orders with respect to T and Vo by standard results in the theory of ordinary differential 
equations (see |l0j|). 

Since r — > oo if and only if the right-hand side of equation (|A2J) vanishes, we define 

T(T, Vo) := 1 - V(T, V ) 2 + ^A(T, V ) 2 (A3) 

and discuss the equation T(T, Vo) = 0. 

Let us observe that T(0, 1) = 0. That is, F = oo for a particle whose initial velocity at 
X = is the speed of light. We will show that there is a unique implicit solution Vo(T) 
(i.e. T(T, V (T)) = 0) such that V (0) = 1 and V (T) < 1 for T near T = 0. Moreover, for 
T > and near T = 0, the implicit solution can be inverted so that T(T(Vo), Vo) = 0; that 
is, for each extremely ultrarelativistic solution, there is a finite time — previously denoted as 
T end — depending on the initial velocity Vo such that r — > oo. Also, T(V ) is a decreasing 
function of Vb such that T(Vq) — > as Vo — > 1. 

The proof of these statements is simply an application of the implicit function theorem. 
Note that T(0, 1) = and Ty o (0, 1) = —2; therefore, by the implicit function theorem, there 
is a function V (T) defined for T near T = such that T(T, V (T)) = and if T(T*, V*) = 
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for (T*, Vq) near (0,1), then V (T*) = V *. To demonstrate that for T > the curve 
(T,Vo(T)) lies in the physical region (Vq(T) < 1), we will show that the Taylor series of Vo 
centered at T = is given by 

V (T) = 1 - -i(GM) 3/2 r " 9/2 T 3 + 0(T 4 ). (A4) 
v2 

In fact, the remaining claims are all simple consequences of this result. 
To establish equation (|A4|) . we first note that 

r (T) = (r 3/2 + ^v / 2GMT) 2/3 (A5) 

and A(0, 1) = 0. A simple computation yields T^(0, 1) = 0, and, from differentiating 
T(T,Vb(T)) = 0, that is, 

T T (T, H(T)) + Ty (T, Vb(T))Vo(T) = 0, (A6) 

we conclude that V^'(0) = 0, where the prime denotes differentiation of the function Vq. 
Proceeding in the same manner to compute the higher-order derivatives of Vq(T), we find 
that Vq(0) = and then differentiating equation (|A6|) twice leads to 

C(0) = ^T TTT (0,1). (A7) 

Next, differentiating equation (|A1|) twice with respect to time leads to 

X TTT (0, 1) = X TTTT (0, 1) = 2W, (A8) 

where W is given by 

Moreover, it follows from differentiating equation (jA3|) three times with respect to T that 

T TT (0, 1) = 0, T TTT (0, 1) = -2X TTTT (0, 1) + 6W. (A10) 

Therefore, X TTTT (0, 1) = T TTT (0, 1) = 2W and V"'(0) = W, as required. 

Starting from equation ()A4|) . we can now derive a simple approximate relationship be- 
tween To and T en( j/(GM) that is valid for extremely ultrarelativistic particles: 

T^y 1 (GMyOft 2 



MM) - V2 V ro I T <>- < A11 > 
Thus for ro = 100 GM and Tq = 300 we find T en <i/{GM) w 20, in agreement with the 
numerical results of Fig. El 
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